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We propose a nanoscale heat engine that utilizes the physics of resonant tunneling in quantum 
dots in order to transfer electrons only at specific energies. The nanoengine converts heat into 
electrical current in a multiterminal geometry which permits one to separate current and heat flows. 
By putting two quantum dots in series with a hot cavity, electrons that enter one lead are forced 
to gain a prescribed energy in order to exit the opposite lead, transporting a single electron charge. 
This condition yields an ideally efficient heat engine. The energy gain is a property of the composite 
system rather than of the individual dots. It is therefore tunable to optimize the power while keeping 
a much larger level spacing for the individual quantum dots. Despite the simplicity of the physical 
model, the optimized rectified current and power is larger than any other candidate nano-engine. 
The ability to scale the power by putting many such engines into a two-dimensional layered structure 
gives a paradigmatic system for harvesting thermal energy at the nanoscale. We demonstrate that 
the high power and efficiency of the layered structure persists even if the quantum dots exhibit some 
randomness. 

PACS numbers: 73.63.-b,85.80.Fi,85.35.-p,84.60.Rb 



I. INTRODUCTION 



Energy harvesting is the process by which energy is 
taken from the environment and transformed to provide 
power for electronics.^ Specifically, thermoelectrics can 
play a crucial role in future developments of alternative 
sources of energy. Unfortunately, present thermoelectric 
engines have low efficiency.^ Therefore, an important task 
in condensed matter physics is to find new ways to har- 
vest ambient thermal energy, particularly at the small- 
est length scales where electronics operate. Utilizing the 
physics of mesosopic electron transport for converting 
heat to electrical power is surprisingly a relatively re- 
cent endeavor. While the general relationships between 
electrical and heat currents and their responses to ap- 
plied voltages and temperature differences have been un- 
derstood since the work of Onsager,^ the investigation 
of thermoelectric properties, and in particular the de- 
sign of nano-engines has taken the form of several fairly 
recent concrete proposals. In 1993, Hicks and Dres- 
selhaus investigated the thermoelectric properties of a 
mesoscopic one-dimensional wire.^ Mahan and Sofo sub- 
sequently showed that the best energy filters are also the 
best thermoelectrics.^ This suggests the use of quantum 
dots with discrete energy levels to investigate thermo- 
dynamic questions.^ The Seebeck effect - the appear- 
ance of a voltage when there is a temperature difference 
across a sample - was investigated for a single quantum 
dot with a resonant level by Nakpathomkun et al.^ Res- 
onant levels were also used as energy filters to make a re- 
lated heat engine from an adiabatically rocked ratchet.^ 
Humphrey et al. note their model can be generalized to 
a static, periodic ratchet, which is a quantum version of 
the model with state-dependent diffusion.^ References 8 
and 9 show that a single resonant level is an ideal heat en- 



gine, and investigate power and efficiency in that system 
which has similarities to the present work. 

Coulomb-blockaded dots can be ideally efficient con- 
verters of heat to work, both in the two-terminaP^' and 
three-terminal^^ cases; however, since transport occurs 
through multiple tunneling processes, the net current 
and power are very small. In light of the small cur- 
rents and power produced by Coulomb-blockaded quan- 
tum dots, open cavities with large transmission that 
weakly changes with incident electron energy have been 
considered.^ While this system produces more rectified 
current than Coulomb-blockaded quantum dots, simply 
increasing the number of quantum channels does not help 
because the energy dependence of transmissions in typi- 
cal mesoscopic conductors is a single-channel effect even 
for a many-channel conductor. Hence, the rectified cur- 
rent and power drop as channel number is increased un- 
less special engineering of the contacts is made. Con- 
sequently, we should optimize the rectified current and 
electrical power for a strongly non-linear system, operat- 
ing in the single channel limit. We do this here with the 
physics of resonant tunneling through barriers connected 
to the hot energy source. We show this achieves the max- 
imum current a single channel can give, which we con- 
jecture is the optimal configuration. Importantly, in our 
model the energy source is separate from the electrical 
circuit, as is required for energy harvesting. Therefore no 
charge is extracted from it and the optimal configuration 
is independent of the parameters of the energy s ource. 
This is different from single quantum dot proposal^^'^'^ 
where the leads are held at different temperatures, and 
the amount of transferred heat depends on the lead chem- 
ical potentials. 

Resonant tunneling is a quantum mechanical effect, 
where constructive interference permits an electron tun- 
neling through two barriers to have unit transmission. 
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Cavity 




Carnot efficient stopping configuration 

FIG. 1. Nanoscale heat engine created from a hot cavity con- 
nected to cold reservoirs via resonant tunnehng quantum dots, 
each containing a single relevant energy level at energy El,r. 
The cavity is kept hot through coupling to an energy reser- 
voir at temperature Tc (not shown) which is considered larger 
than the reservoir temperature Tr. Thermal broadening of 
the Fermi functions in the three regions (source, cavity and 
drain) is shown by the light shading, (a) Rectification config- 
uration. In the absence of bias (short circuit), electrons enter 
the cavity via the left lead, gain energy = Er—El from 
the cavity, and exit through the right lead, transferring an 
electrical charge e through the system, (b) Carnot efficient 
stopping (open circuit) configuration. 

This is only true if the electron has a particular energy 
equal to the bound state in the quantum dot, or within a 
range of surrounding energies, whose width is the inverse 
lifetime of the resonant stateP^ Electrons with any other 
energy are effectively forbidden from transmitting across 
the quantum dot. In this way, a resonant tunneling bar- 
rier acts like an energy filter: only those electrons that 
match the resonant condition are permitted to pass. We 
assume for simplicity the resonant tunnel barrier (or the 
dot) is symmetrically coupled. 

The geometry we consider is related to one that has 
already been fabricated experimentally, but considered 
for a different pur pose, based on the theoretical proposal 
of Edwards et al.^^^ Prance et al. performed exper- 
iments on this system in its dual role as an electronic 
refrigerator.'^ They demonstrated that applying bias to 
the system results in cooling a large 6 fim^ cavity from 
280 mK to below 190 mK. We are primarily interested in 
cavities made as small as possible while still having good 
thermal contact with the heat source. This further minia- 
turization permits many engines to be put in parallel and 
give a large output power. 

The paper is organized as follows. In Sec. [lljwe de- 
scribe the transport through our device, the results of 
analytical and numerical calculations are presented in 



Sec. Ill Section IV investigates the scaling of the sim- 



ple system consisting of two quantum dots up to a lay- 
ered structure with many channels contributing at the 
same energies and discusses the robustness of the pro- 
posed structure to fluctuations in the fabrication process. 
The conclusions are presented in Sec. |Vl 



II. MODEL 

The model we consider (shown in Fig. 1) consists of a 
cavity connected to two quantum dots, each with a reso- 
nant level of width 7 and energy El r. We consider the 
situation where the widths are equal, while the energy 
levels are different (these are controlled by gate voltages) . 
The energy difference = Er — El is an important 
energy scale of our composite system, which we refer to 
as the energy gain. It is distinct from the level spacing 
5 in the individual dots as well as from the level width. 
The nano-cavity the dots are connected to is considered 
to be in equilibrium with a heat reservoir of temperature 
Tc that is hotter than the left and right electron reser- 
voirs, having chemical potentials (jll.r and equal tem- 
peratures, Tr. We assume strong electron-electron and 
electron-phonon interactions relax the electron energies 
as they enter and leave the cavity, so the cavity's occu- 
pation function may be described with a Fermi function, 
f{E — jii^T) = l/(l + exp[(T^ — /i)//cBT]) completely char- 
acterized by a cavity chemical potential /ic and temper- 
ature Tc, with /cb being the Boltzmann constant. This 
process of inelastic energy mixing is assumed to occur 
on a faster time scale than the dwell time of an electron 
in the cavity. Thermal energy flows from the coupled 
hot bath into the cavity as a heat current, and keeps 
the temperature different from that of the electron reser- 
voirs. The nature of the heat reservoir is not specified 
in this model, but refers quite generically to any heat 
source we wish to harvest energy from. Our setup should 
be contrasted with the more widely studied two-terminal 
configurations in which electric and heat currents flow 
in parallel. Our model permits a separate heat circuit to 
the central reservoir and a separate (transverse) electrical 
circuit. 

The chemical potential of the cavity and its temper- 
ature (or equivalently, the incoming heat current) are 
constrained by conservation of global charge and en- 
ergy. These constraints are given by the simple equa- 
tions, /L+^i?=0, and Jl-\-Jr-\-J=0 in the steady state, 
where Il,r^^ electrical current in the left or right contact, 
and Jl,r the energy current. Energy current is seemingly 
not conserved because of the heat current J flowing from 
the hot reservoir. 

The currents /j, j=L^R^ are given by the well 
known formulas Ij = {2e/h) J dE Tj{E)[fj—fc] and Jj = 
{2/h) f dETj{E)E[fj- fc], where Tj{E) is the transmis- 
sion function of each contact for each incident electron 
energy E. In our quantum dot geometry, the resonant 
levels give rise to a transmission function of Lorentzian 
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shape 



T1T2 



(1) 



where ri^2 are the attempt frequencies of the two barriers 
of the resonant quantum dot(s). Here we assume sym- 
metric couphng for simphcity, Fi = r2 = 7, so 7 is the 
width of the level, or inverse lifetime of an electron in the 
dot. Note that the Lorentzian energy dependence applies 
if the level width 7 is small compared to the level spacing 
d in the individual quantum dots. A crucial advantage of 
our setup which we will exploit later on is that it permits 
the use of small dots with a level spacing that is large 
compared to temperature but with the energy gain and 
level width of each dot on the order of the temperature. 

In the limit where the width of the level is smaller 
than the thermal energy in the cavity /dot system, 
7<CA:bTc, /c^Tr, the transmission will pick out only the 
energies or in the above energy integral expres- 
sions for the currents giving simple equations. Conse- 
quently, we have these equations for the conservation 
laws for charge and energy: 

= /l - /CL + /i? - fcR. (2a) 
= JV(27) + ElIJl - fcL] + ErI/r - fcRi (2b) 

where = J^El-^l.T^), /r = f{ER-^R,T^), fcL = 
/(£^L-Mc,^c), and fcR = /(£^i?-/ic, ^c)- From these 
two equations, we can solve for (say) the quantity fcR — 
/r = Jh/{2'yAE). This quantity is proportional to the 
electrical current through the left lead II = —Ir = /, 
the net current flowing through the system. 

A solution of Eqs. (2a) and (2b) to linear order in the 



deviation of the cavity's temperature and chemical po- 
tential from the electronic reservoirs indicates that the 
maximal power of the heat engine will be produced when 
the chemical potentials of the reservoirs are symmetri- 
cally placed in relation to the average of the resonant 
levels, /ii^^L=±M/2+(^L+^i?)/2. For this special case, 
an exact solution is possible because the constant solu- 
tion iic = {El^Er) /2 for the cavity chemical potential 
satisfies the charge conservation condition (2a) for all 
temperatures. 



III. RESULTS 
A. Limit of small level width 

We now describe the physics of this nano-engine. We 
first focus on the regime 7 <C A^b^r^^b^c, which can 
be analyzed analytically and will afterwards discuss the 
regime 7~/cbTr, k^Tc which we numerically find to yield 
the largest current and power. Physically, if an electron 
comes in the left lead at energy El and exits the right 
lead with energy Er>El^ it must gain precisely that 
energy difference /S.E = Er — El. Thus, in the steady 
state, any incoming heat current J must be associated 
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FIG. 2. (a) Scaled maximum power as a function of energy 
gain /\E for AT = T and level width 7 and iir optimized 
to give maximum power, (b) Scaled maximum power as a 
function of 1ST jT for optimized values of Er — /\E/2^'^ 
and iiR. (c) Efficiency at maximum power for the val- 
ues of Er = A^/2,7 and fiR chosen to maximize power, 
(d) The optimized values are plotted versus AT/T. Here, 
Er = AE/2, fiR = /i/2. 

with an electrical current /, with a conversion factor of 
the energy gain, AE^ to the quantum of charge, e. 

This results holds regardless of what bias is applied or 
what the temperature is. 

The efficiency of our heat engine, 7^, is defined as the 
ratio of the harvested electrical power P=\{jiiL—I^R)I\/e 
to the heat current from the hot reservoir, J. For our 
system it takes a particularly simple form. 



AE 



(4) 



In order to proceed further, we must find the chemi- 
cal potential of the cavity and its temperature given in 
terms of the incoming heat current and chemical poten- 
tials and temperature of the electron reservoirs. These 
are found by employing the principle of conservation of 
global charge and energy, see Eqs. (2a) and (2b). 



Rather than considering the cavity temperature Tq as 
a function of the heat current J we are harvesting, we 
can turn our perspective around, and keep the cavity 
temperature Tq fixed from being in thermal equilibrium 
with the hot energy source. With this insight, we can 
express the heat current J in terms of the hot cavity 
temperature and other system parameters. We find in 
the limit where 7 <C Aj^Tr, ksTcAE, 

2-fAE^ 



J 



h 



-[/(AE/2,Tc) - f{AE/2 - /i/2,TR)], (5) 



which satisfies charge and energy current conservation. 
In Eq. ([5|, is Planck's constant. 

Importantly, without bias there is a rectified electrical 
current given by 

I = eJ/AE ^ ^-l^[(kBT^)-' - {kBTc)-\ (6) 
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in the limit where kBTji^kBTc ^ AE. This current is 
driven solely by the fixed temperature difference between 
the systems. We note both the heat and electrical current 
are proportional to 7, the energy width of the resonant 
level. Consequently, the currents and power produced 
in this system will tend to be small since we have as- 
sumed that 7 is the smallest energy scale. It is also clear 
that both are controlled by the size of AE^ so increasing 
this energy gain will improve power until it exceeds the 
temperature. Later we will generalize these results by 
numerically optimizing the power produced in this nano- 
engine. 

In order to harvest power from this rectifier, a load 
should be placed across it. Equivalently, we could apply 
a bias V = /i/e to this system tending to reduce the rec- 
tified current. At a particular value, /istop, the rectified 
current vanishes, giving the maximum load or voltage one 
could apply to extract electrical power at fixed tempera- 
tures Tr,Tc. This value is found when J and / vanish, 
given by Eq.([5|: 



Mstop 



AE 1 



Tc 



(7) 



Consequently, the voltage applied must not be larger 
than /istop/e and therefore from Eq. Q, the efficiency 
is bounded by 7^ < 1 — ^ = 77c. At the stopping volt- 
age, the thermodynamic efficiency attains its theoretical 
maximum, the Carnot efficiency, r]c , showing this system 
is an ideal nanoscale heat engine. Naturally, at this point 
[see Fig. [ijb)] the system is reversible with no entropy 
production. Also interesting is the efficiency at the bias 
point where power is maximum. For temperature larger 
than or eV=ii^ we can approximate the Fermi func- 
tions to find P ^ (7/4/i/cBTR)/i(/istop— /i), resulting in a 
parabola as a function of /i, with maximum power 



(8) 



and efficiency T^maxP = which is in agreement with 
general thermodynamic bounds for systems with time- 
reversal symmetry.l^M^ 



B. Optimization 

One can go beyond this limit for the efficiency by solv- 
ing the conservation laws numerically. We optimize the 
total power produced by the heat engine by varying the 
resonance width 7, as well as the energy gain AE and 
applied bias F=/i/e, given fixed temperatures Tr,Tc. 
These results are plotted in Fig. [2j where we define the 
average temperature T=(Tc+Tr)/2, and its difference, 
AT=Tc-Tr. In Fig. 2(a), we see that for AE < ksT 
the power increases as A£^^, as indicated in Eq. ([8|, 
but then levels off and decays exponentially, attaining 
its maximum around AE = ^ksT. Similarly, the choice 
7 = ksT gives optimal power. We emphasize that AE 
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FIG. 3. Scaled heat current J/[2{kBTf/h] leaving (blue) or 
entering (red) the central cavity vs temperature difference (x 
axis) and applied bias (y axis) divided by average tempera- 
ture. The plot is for system parameters optimized for maxi- 
mal power output, an energy gain of AE ^ GksT and level 
width 7 ^ ksT. The green line is the J = curve for 7^0 
while the black curve is the J — line for the optimized 7. 
The system can work as a heat engine (HE) in the blue re- 
gion for AT > — 11/2 > (configuration shown in the 
inset labeled HE), or as a refrigerator (R) in the blue region 
for AT < = /i/2 < (configuration shown in the inset 
labeled R). 

and the level width are two essentially independent en- 
ergy scales in our system. As a matter of fact the energy 
gain AE ^ ^ksT is almost an order of magnitude larger 
than the level width 7 ^ ksT. These considerations sug- 
gest an experimental strategy for maximizing the power 
of such a device: Measure what the resonant level widths 
are, and tune the reservoir temperatures and energy gain 
(with the help of gate voltages), to it. 

From Fig. ^h) we also see that the efficiency at max- 
imum power drops from half the Carnot efficiency to 
about 0.27^c when the parameters are optimized. How- 
ever, we note that when 7 is kept small, in the nonlinear 
regime the efficiency can exceed the bound r/maxP ^ 
found in the linear regime. This small drop in efficiency 
is more than compensated by the extra power we ob- 
tain. Importantly, according to Fig. |2j the power reaches 
a maximum of Pmax'^0-4(/cBAT)^//i, or about 0.1 pW 
at AT=1K, a two order of magnitude increase from 
a weakly nonlinear cavity. This jump in power can 
be attributed to the highly efficient conversion of ther- 
mal energy into electrical energy by optimizing both the 
level width and energy level difference. Compared to a 
heat engine based on resonant tunneling through a sin- 
gle quantum dot in a strictly two-terminal geometry,!^ 
our three-terminal energy harvester gains a factor of 2 
in power while achieving the same efficiency. To give 
some perspective on this output, if a 1 cm^ square array 
of these nanoengines were fabricated, each occupying an 
area of 100 nm^, they would produce a power of 0.1 W, 
operating at AT = 1 K. 

In Fig. [Sj the heat current J is plotted versus tem- 
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FIG. 4. Self-assembled dot engine: A cold bottom electrode 
(dark blue) is covered with a layer of quantum dots (orange) 
embedded in an insulating matrix (transparent blue). The 
quantum dot layer is covered by the hot central region shown 
in dark red. On top of it, there is another quantum dot layer. 
The whole structure is terminated by a cold top electrode 
(dark blue). Importantly, the positions of the quantum dots 
in the two layers do not have to match each other. Thus, 
the device can be realized using self- assembled quantum dots. 
The latter can have charging energies and single-particle level 
spacings of the order of lOmeV (Refs. 21 and 22), thereby 
allowing the nanoengine to operate at room temperature. 

perature difference and applied bias. There we find that 
when system parameters are optimized to give maximum 
power, the system can be operated in the mode of a heat 
engine (HE) or a refrigerator (R). However, in contrast 
to the case where the levels are narrow compared to the 
other energy scales [and consequently the cavity can cool 
to arbitrarily low temperatures in principle; see green 
solid line, Eq. ([7|)], for this choice of parameters the cav- 
ity will only cool to the temperature where the J = 
curve (solid black line) bends back. 



C. Scaling 

In reality, there will be other quantum dot resonant 
levels the electron can occupy that are higher up in en- 
ergy. We have assumed that the cavity temperature and 
applied bias are sufficiently small so transport through 
these levels can be neglected. Our model is quite general, 
so it may be applied to both semiconductor dots in two- 
dimensional electron g ases,^*^ as well as three dimensional 
metallic dots.f^ This latter case is quite interesting since 
one can fabricate an entire plane of repeated nanoengines 
in parallel in order to scale the power. Furthering this 
idea, for such a repeated array of cavities and quantum 
dots, one can connect all the cavities to make a single 
engine, see Fig. [4j The two boundary layers consist of 
planes of quantum dots, so electrons can only penetrate 
through them. These layers sandwich a hot interior re- 
gion and separate it from the left and right cold exterior 
contacts. This is equivalent to taking a large cavity with 
two leads and scaling the power by adding more quantum 



dots (rather than trying to add more channels to a single 
contact). Interestingly, the layered structure can help to 
reduce phononic leakage heat currents that would other- 
wise reduce the efficiency of the system: Phonons scatter 
efficiently at interfaces, and the random dot arrangement 
will further reduce phonon related heat losses. The sand- 
wich engine fabricated with self-assembled quantum dots 
tolerates variations in width and fluctuations in energy 
levels, as we will see in the next section. 



IV. ROBUSTNESS OF THE LAYERED 
QUANTUM DOT ENGINE TO FABRICATION 
FLUCTUATIONS 

Here, we further investigate the self assembled quan- 
tum dot layered engine, and what some of its theoretical 
characteristics are under realistic fabrication conditions. 
The basic operating configuration for the engine is shown 
in Fig. |5j Heat flows from the hot energy source we wish 
to harvest energy from into the engine, and is converted 
into electrical power, along with residual thermal energy, 
dumped into the cold temperature bath. The electri- 
cal current is carried by the cold thermal bath, where it 
powers a load and then completes the electrical circuit on 
the opposite cold terminal. The flow of heat out of the 
hot energy source will consequently tend to cool the hot 
source. The energy harvesting application we primarily 
have in mind is taking heat away from computer chips 
and running other devices on the chip itself. In electrical 
chips, thermal energy is an abundant and free resource. 
Indeed, heat is not only free, it is a nuisance preventing 
further improvements on chip technology. The fact that 
the proposed heat engine not only harvests the thermal 
energy, converting some of it into electrical power, but 
also cools the hot source is therefore a side benefit to the 
proposal of including nanoscale heat engines as part of 
an emerging chip technology. 

In the actual fabrication of a resonant tunneling nano- 
engine, as long as there are only a few dots, the precise 
placement of the resonant levels can be controlled by gate 
voltages in order to maximize the power generated by the 
engine. However, as soon as we consider self- assembled 
quantum dots with charging energies and single-particle 
level spacings of the order of 10 meV,l^^^^ thereby allow- 
ing the nano-engine to operate at room temperature, this 
kind of control is out of the question. To make such an 
engine, there are several possible fabrication techniques 
that could be employed using layers of quantum dots and 
wells to have the resonant energy levels lower than the 
Fermi energy on one side of the heat source, and higher 
on the other side. However, in all of these fabrication 
methods, the growth of quantum dots does not occur at 
a perfectly regular rate, so it is natural to expect there 
will be variation of the resonant energy level from dot to 
dot. We must then check whether this fact will degrade 
the performance of the engine, and if so by how much. 

To answer this question, we consider the energy- 
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FIG. 5. Schematic of the system operation configuration. 
Heat enters the engine from the pink hot region, indicated 
by the red chevrons. The engine itself is signified by the yel- 
low sandwich with black holes indicating the position of the 
resonant tunneling quantum dots. The position of the holes 
can be disordered, similar to a slice of Swiss cheese. Electri- 
cal current is generated perpendicular to the layers, indicated 
with the blue arrows flowing in the light blue cold region. 

resolved current arising from a number of electrons pass- 
ing through N quantum dots on the left and then the 
right layer. The total current coming from the left slice 
is given by 



rtot 

r. 



2e 

N 



j dET,^^L[f{E - fiL.Tn) - f{E - ^c.Tc)]^ 



(9) 



with similar equations for the total right current, as well 
as the energy currents. Here, Ti{E^Ei) is the transmis- 
sion probability of quantum dot i, which has a resonant 
energy level Ei and a width 7^, Ti{E,Ei) = jf/[{E — 
Ei)'^-\-jf] for symmetrically coupled quantum dots. Since 
neither the left nor cavity Fermi functions depend on the 
level placement, the sum over the quantum dots can be 
done to give an effective transmission function Tq^^^l for 
the whole left slice. We can make further progress by 
assuming the fabrication process can be described as a 
Gaussian random one, where the energy level Ei is a 
random variable with an average of E^ and a standard 
deviation of a. For simplicity, we only consider random 
variation in Ei^ but there will also be variation in 7^ we 
ignore for the present. With this model, the effective 
transmission will have the average value 



(Teff.i) = N{Ti{E,E))p = N / dET{E 



/- 



E)PG{E,a), 
(10) 
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FIG. 6. The power per nanoengine is plotted vs the width of 
Gaussian random distribution of energy levels. The power is 
normalized to the maximum power a single nanoengine can 
give (for optimized parameters), while the distribution width 
is plotted in units of A^, the average energy gain between 
left and right dots, (a) The power is plotted for different level 
widths of the dots, where 7opt = 1.02 ksT. (b) The power is 
plotted for different applied voltages, where /in, opt = 0.7 ksT. 



where Pq is the Gaussian distribution described above. 
Thus, we see the effective transmission function is sim- 
ply a convolution of the Lorentzian transmission function 
and the Gaussian distribution, known as a Voigt pro- 
file. This leads to further broadening of the Lorentzian 
width. With these considerations, the conservation laws 
for charge and energy retain the same basic form as de- 
scribed in Eqs. (2a) and (2b), but with N times the Voigt 



profile playing the role of the energy-dependent trans- 
mission for the left and right leads. We have numerically 
solved these equations and plotted the maximum power 
per nanoengine versus the width of the Gaussian distri- 
bution in Fig. [6j The parameters are chosen so as to 
optimize the engine's performance without any random- 
ness in the level position. As the randomness of the level 
position is increased, the power begins to drop as ex- 
pected. However, even when the scatter of the energy 
levels is 10% of AE", the power only drops to 90% of its 
maximum, showing that this engine is robust to these 
kinds of fluctuations in fabrication. Notice that if the 
level width is less than the optimal amount, some disor- 
der in the level energy can actually improve performance 
in comparison to a cleanly optimized level width. This 
is because of the additional broadening in energy space 
the level disorder provides. Another interesting effect is 
shown in Fig. [7| where the value of 7 is not optimal to- 
gether with a given amount of disorder in the energy level 
positions. The figure demonstrates that even in this ex- 
perimentally realistic case, a change of the voltage and 
relative level spacing can yield results which are nearly 
as good as the optimal case. 
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FIG. 7. The power per nanoengine is plotted vs the right 
chemical potential and the average energy gain, A^, in units 
of the optimal power, Pmax- The level width is taken to be 
a suboptimal value, 7 = O-SksT and the width of Gaussian 
random distribution is fixed at a = 0.2kBT. The subscript 
"opt" refers to the parameter that is optimized for maximal 
power in the clean case. 



V. CONCLUSIONS 



In this paper, we have opened a route to highly effi- 
cient solid state energy harvesting. Our work shows that 
the most efficient structure is found with energy filters 
that have transmission probabilities close to 1 in a con- 
figuration where the level structure of spatially separated 
quantum dots can be independently adjusted to give the 
desired energy gain. This yields a high power heat engine 
with no moving parts. We show that nano-scale quan- 
tum dots can be employed in a parallel configuration that 
delivers substantial power. 
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